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Correct the first paragraph on p 5 to read: 

Instead of referring the reader to earlier papers on this subject 
(e.g., ref. 1 +) , a brief presentation of the procedure used to calculate 
the detonation parameters of carbon monoxide - air mixtures is given 
here. 

The correct list of references (p 14 ) is as follows: 
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SUMMARY 


/O / 

The composition, temperature, pressure and density behind a stable 
detonation wave and its propagation rate have been calculated for five 
mixtures of carbon monoxide and air at 1, 5, 25 and 100 atmospheres 
initial pressure, and at an initial temperature of 313 *l 6 °K. The results 
of these calculations show that the detonation velocities of stoichio- 
metric carbon monoxide - air mixtures increase with increasing ini- 
tial pressure whereas those of very rich and very lean mixtures are 
independent of the initial pressure. 

The pressure ratios across the wave and the temperatures behind the 
wave increase with increasing pressure for near -stoichiometric mixtures, 
but they are unaffected by pressure for rather rich and lean mixtures. 

INTRODUCTION 

Complete thermodynamic and chemical equilibrium is assumed to exist 
behind the detonation wave. Dissipating effects such as heat transfer, 
viscosity, diffusion, and chemical reaction rates have been disregarded 
since it was not intended to determine the structure of the wave. Details 
of the method used for calculating the detonation parameters were pub- 
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lished previously (ref. l), and calculations for various hydrogen-oxygen 
mixtures over the same range of initial pressures can be found in 
reference 2. 

The detonation parameters have been derived from the Hugoniot 
equation for the reacted gas mixture in equilibrium, and from the Chap- 
man -Jouguet condition that the detonation velocity is the minimum wave 
velocity. Therefore, the results are based on the Chapman -Jouguet 
point at which the velocity of the reacted gas (relative to the detona- 
tion wave) is equal to the equilibrium sonic speed (ref. 3)« 

Besides being of fundamental interest, the theoretical detonation 
parameters, particularly the detonation velocity, is of importance in 
the study of the formation of detonation waves. These values are used 
to determine the detonation induction distance. Results from these 
experiments will be published separately. 


SYMBOIS 


% 

AE 

o 

h a 



sonic velocity in gas entering shock wave 
heat of reaction at absolute zero temperature 
absolute enthalpy of gas mixture per unit mass leaving 
shock wave 

absolute enthalpy of gas mixture per unit mass entering 
shock wave 

heat of formation of species i at temperature T & formed 
from elements, making species i, in the standard state 
(p = 1 atm, T = 298.16 °K) 

heat of formation of species i at temperature formed 
from elements, making species i, in the standard state 
(p = 1 atm, T = 298.16 °K) 
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dimensionless heat of formation of species i at temperature 
T formed from elements , making species i, at standard 
state (p = 1 atm, T = 298. l6 °K) 

dimensionless heat of formation of species i at T a formed 
from the elements, making species i, at standard state 
(p = 1 atm, T = 298.16 °K) 

dimensionless heat of formation of species i at formed 
from elements, making species i, at standard state 
(p = 1 atm, T = 298.16 °K) 

equilibrium constant of reaction j based on partial pressures 

Mach number of gas mixture entering shock vave 

Mach number of stable detonation vave . Mp. = — 

% 

mole number of gas mixture, per unit mass of mixture, entering 
shock vave 


n 


mole number of species i, per unit mass of mixture, entering 
shock vave 


K total mole number of element k in unit mass of mixture 

k 

p absolute pressure of gas behind shock vave 

a 

p absolute pressure of gas behind stable detonation vave 

a,U 

p _ _ absolute pressure of gas behind normal shock vave vithout 

& , I. . s . 

chemical reactions 
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p. 

i,a 


P imp 



(ft 

T 


T 

a,N.S. 


AT 


U 


U, 


14. 


u 


a,D 


v 


7b 


P 


a,D 


P 


b 


absolute partial pressure of species i behind shock wave 
absolute pressure of gas entering shock wave 
absolute impact pressure of detonation vave 


total quantity of element k, in atomic state, in reacted 

r: k ir 


mixture in terms of pressure units, = N, 


universal gas constant 

absolute temperature of gas behind shock wave 
absolute temperature of gas behind detonation wave 
absolute temperature of gas behind normal shock wave 
without chemical reactions 
absolute temperature of gas entering shock wave 
temperature difference between estimated temperatures 
linear velocity of gas leaving detonation or shock wave 
linear velocity of gas entering shock wave 

linear velocity of detonation wave relative to unreacted gas 

linear velocity of detonation wave relative to reacted gas 

specific volume, v = i 

P 

isentropic coefficient (ratio of specific heats) of unreacted 
gas mixture 

density of gas mixture behind detonation wave 
density of gas mixture entering shock wave 
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METHOD OF CALCULATION 


Instead of referring the reader to earlier papers on this subject, 
a brief presentation of the procedure used to calculate the detonation 
parameters of carbon monoxide - air mixtures is given here. 

It has been assumed that only COg, CO, O 2 , 0, N 2 and NO occur to 

an appreciable extent in the detonation wave of carbon monoxide - air 
(0.79 mole Ng plus 0.21 mole 0^) mixtures ranging in fuel concentration 

from 10 per cent CO to 50 per cent CO. The presence of atomic nitrogen, 
higher oxides of nitrogen, free carbon, and electrically charged species 
has been ignored because of the extremely low concentrations found for 
these species in the mixtures considered during the present study. Thus 
the composition of the detonated gas is governed by the following six 
equations : 


(C0 2 i=* CO + | 0 2 j = C0 2 ) (1) 
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(6) 


g V P NO = IP M 

p„^ + P„^ Ip N„ 

^C0 2 *C0 L c 


p co 2 + P C0 + 


'Or 


+ + P„. + P„„ = 2 P 4 _ = P 


’N 2 


NO 


i,a 


For a temperature, T , in the vicinity of the expected detonation 

81 

temperature and for an arbitrarily assumed value of p the partial 

°2 

pressure of molecular nitrogen can be calculated by using equations (l) 
through (5). 



(T) 
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The remaining partial pressures follow easily. From equation (2) 


P« = . From equation (3) p = *\ /~p~ «~p ’ • K . From 

o y o 2 o 2 NO \J o 2 n 2 NO 

equations ( 5 ) and (l) 



And from equation (l) 



For these partial pressures and with the first estimate of the detona- 
tion temperature we obtain from the Hugoniot equation 
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(8b) 


and from equation (6) we also obtain a value for this pressure ratio 


r t> p +p +p +p+p +p 

P i,a y C0 2 P C0 P 0 2 P 0 N 2 P N0 


(6) 


If the pressure ratios derived from equation (8b) and equation (6) 
differ, the calculations must be repeated with a new value of p 

U 2 

for the same value of T , until agreement is obtained. These 

a,l 
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calculations are carried out for three values of T (T , < T ^ < T ). 

a v a,l a , 2 a, 3' 

The three corresponding Mach numbers (multiplied by the specific heat 
ratio 7 ) are calculated from the relationship 






E P, 


l2 a 


1 



^ n> % 


( 9 ) 


The Mach number of the stable detonation wave is obtained by locating 

O 

the minimum of the parabola which can be drawn when the three 7^ 

values are plotted versus the three assumed T a values. To obtain 

maximum accuracy with this interpolation method, the temperatures 

T -t , T 0 and T 0 must be selected in such a fashion that 
a, l 7 a, 2 a, j 

(io) 


If equal temperature intervals between the three assumed temperatures 
are chosen so that 




a,l 



T a,l } 


AT 


( 11 ) 
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the parabolic interpolation method leads to 
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Then, according to equation (9), the detonation velocity becomes 


U D 


-i 


7 M 
b D 


n 


& T > 


(13) 


For the temperature behind the detonation wave, the interpolation 
yields 


T 


a,D 



AT 

2 


4T [<\ \\ - Oj 

[< 7 b Ol - < 7 b %%] + [ (7 b M b 2) 3 - (7 b V%] 


(14) 


Finally, values of p _ and p o n are obtained by linear interpolation 

&,) JJ 

between the respective values calculated for the three temperatures 
chosen. If the temperature intervals, AT, are sufficiently small 
(e.g., 10°K), the error due to this linear interpolation is negli- 
gible . 

The impact pressures are calculated according to the relation- 
ship 


V, 


imp 


= p - p + p (U - U J 2 = ( p 


a,D 


a,D 


D 


a,D 


pj (15) 
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The pressure and temperature behind the normal shock which forms 
the front of the detonation wave have been calculated on the basis that 
no chemical reactions occur in the shock zone and that the vibrational 
degrees of freedom in the molecules of the initial gas are inactive. 

Thus the usual equations employing constant specific heats could be used 
for calculating the conditions behind the normal shock propagating at 
the same rate as the detonation velocity. 


DISCUSSION OF RESULTS 

The results of the calculations Eire compiled in tables 1 through 
5 and they are depicted graphically in figures 1 through 6. Values 
of the equilibrium constants used for the calculations are given in 
table 6 and the dimensionless heats of formation are listed in table 
7. Both thermodynamic functions have been tabulated only for large 
temperature intervals to save space. These data represent carefully 
checked and, in some instances, recalculated values of the statistical 
calculations of several authors (refs. 5-14). A method to correct old 
values with respect to new data on atomic and molecular constants has 
been described in reference 15 . Values for 10°K intervals were obtained 
by interpolation. 

From figure 1 it can be seen that the detonation velocity increases 
only slightly when the initial gas pressure of the combustible mixture 
is raised. For example, consider a carbon monoxide - air mixture 
containing 29*6 per cent carbon monoxide. For this mixture the deto- 
nation velocity increases only from 1668 m/sec to 1744 m/sec for a 
rise in the initial pressure from 1 to 100 atmospheres. The maximum 
of the curves depicting the velocity as a function of the initial gas 
composition occurs in the region of rich mixtures at low initial 
pressures. As the initial pressure is increased the maxima of these 
curves approach the stoichiometric fuel concentration because at the 
higher pressures less dissociation occurs in the reacted gas of the 
wave. In the case of carbon monoxide - air mixtures the sonic speed 
in the initial gas is practically independent of its composition. 
However, in hydro gen -oxygen mixtures the sonic speed increases rather 
markedly when the concentration of hydrogen is increased. Therefore, 
the maximum detonation velocity of hydrogen-oxygen mixtures (ref. 2) 
occurs in rather rich mixtures even at very high pressures. The 
small reduction of the degree of dissociation in this range of fuel 
concentrations is not as effective as the increase in the speed of 
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sound in the initial gas mixture. As a result of this condition the 
maximum of the curves representing the detonation velocity as a 
function of the hydrogen concentration for hydro gen -oxygen mixtures 
shifts only very slightly towards the stoichiometric fuel concen- 
tration as the initial gas pressure is increased (see fig. 3> ref. 2). 
That the detonation velocity of a combustible gas mixture is propor- 
tional to the speed of sound in the unburned gas is indicated by 
equation (13). 

Figure 2 illustrates the effect of initial pressure on the 
pressure ratio across the detonation wave. For a mixture containing 
29.6 per cent carbon monoxide, this ratio increases by about 10 per 
cent when the pressure of the initial gas is increased from 1 to 100 
atmospheres. Following the behavior of the detonation velocity, the 
pressure ratio as a function of the fuel concentration also has a 
maximum in the region where the mixtures are rich in fuel and, with 
increasing initial pressure, this maximum also moves toward the 
stoichiometric fuel concentration. According to these calculations, 
the pressure behind the detonation wave of carbon monoxide - air 
mixtures is only 20 per cent lower than that of hydro gen -oxygen 
mixtures having the same equivalence ratio. The pressure ratios 
across the normal shock preceding the detonation wave, figure 6, 
are approximately 80 per cent greater than the detonation pressure 
ratios . 


The curves representing the gas temperature behind the detonation, 
figure 3 ) are quite similar in shape to those depicting the detonation 
velocity, figure 1. When the initial pressure is increased from 1 to 
100 atmospheres, the detonation wave temperature of a mixture containing 
29.6 per cent carbon monoxide increases from 2864°K to 31T3°K. This 
temperature increment amounts to almost 11 per cent of the lower 
temperature. Again, at the lower initial pressures the maximum temper- 
ature is produced by rather fuel-rich mixtures and this maximum 
approaches the stoichiometric fuel concentration as the initial 
pressure is raised. The temperature behind the normal shock preceding 
the detonation wave, figure 5> is approximately ^5 per cent lower than 
the temperature in the fully detonated gas. 

The mole fractions of the various species occurring in the deto- 
nation wave are graphically represented in figures 7 through 10 for 
each of the four initial pressures of 1, 5> 25 and 100 atmospheres. 
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These curves demonstrate clearly the well-known fact that the degree 
of dissociation is less at the higher pressures whereas the absolute 
concentration of free radicals increases with increasing pressure. 

As anticipated, the mole fraction of molecular oxygen decreases with 
increasing carbon monoxide content in the mixtures. However, the 
mole fraction of atomic oxygen attains a maximum in nearly stoichi- 
ometric mixtures. 

The precision of the calculations is about equal to the accuracy 
of the thermodynamic data used. Because of the somewhat large temper- 
ature intervals used in the present calculations, the estimated 
precision attained is about 0.1 per cent of the value of the calcu- 
lated parameters . 


CONCLUSIONS 

Although the present method of calculating detonation parameters 
is intended primarily for use with desk calculators, it is too slow 
when very many calculations are to be made. However, the method can 
be adapted for electronic computing machines quite easily. 

The detonation parameters of carbon monoxide - air mixtures are 
affected less by changes in the initial pressure than those of 
corresponding hydro gen -oxygen mixtures because of the lower degree 
of dissociation of the former. 


The Ohio State University 

Rocket Research Laboratory 

Columbus 10, Ohio May 15, 1962 


This investigation was conducted under Grant NsG -44-60 at 
The Ohio State University under the sponsorship and with the 
financial assistance of the National. Aeronautics and Space 
A dmi ni strati on . 
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TABLE 1. COMPOSITION OF COMBUSTION GAS IN DETONATION WAVE* 


Mo lax 
per cent 
fuel in 
mixture 

P b 

atm 

p co 2 

P a,D 

p co 

P a,D 

p o 2 

P a,D 

p o 

P a,D 

P a,D 

P 

NO 

p a,D 

10 

1 

0.1051 

0.0000 

0.1456 

0.0000 

0.7475 

0.0018 


5 

0.1051 

0.0000 

0.1456 

0.0000 

0.7475 

0.0018 


25 

0.1051 

0.0000 

0.1456 

0.0000 

0.7475 

0.0018 


100 

0.1051 

0.0000 

0.1456 

0.0000 

0.7475 

0.0018 

20 

1 

0.2121 

0.0090 

0.0724 

0.0013 

O.69I8 

0.0134 


5 

0.2165 

0.0050 

0.0706 

0.0007 

0.6930 

o.oi4i 


25 

0.2194 

0.0025 

O.0694 

0.0003 

O.6938 

0.0145 


100 

0.2207 

0.0013 

O.O689 

0.0002 

0.6943 

0.0146 

29.6 

1 

0.2578 

0.0757 

0.0294 

0.0033 

0.6202 

0.0135 


5 

0.2735 

0.0625 

0.0232 

0.0020 

0.6247 

0.0140 


25 

O.2890 

0.0495 

0.0173 

0.0011 

0.6294 

0.0138 


100 

0.3011 

0.0392 

0.0128 

0.0006 

0.6333 

0.0130 

bo 

1 

O.2615 

0.1904 

0.0073 

0.0020 

0.5322 

0.0067 


5 

0.2737 

0.1809 

0.0034 

0.0008 

0.5361 

0.0051 


25 

O.2817 

0.1746 

0.0011 

0.0003 

0.5391 

0.0032 


100 

0.2853 

0.1718 

0.0004 

0.0001 

0.5407 

0.0019 

50 

1 

0.2324 

0.3257 

0.0004 

0.0002 

0.4403 

0.0011 


5 

0.2339 

0.3246 

0.0001 

0.0001 

0.4409 

0.0005 


25 

0.2344 

0.3242 

0.0000 

0.0000 

0.4411 

0.0002 


100 

0.2346 

0.3240 

0.0000 

0.0000 

0.4412 

0.0001 


* Initial temperature = 313.l6°K 
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TABLE 2. PROPERTIES OF A CARBON MONOXIDE - AIR DETONATION WAVE, p = 1 atm 

b 


Molar 
per cent 
fuel in 
mixture 

T b 

°K 

T a,D 

°K 

T a,N.S. 

°K 

p a,D 

atm 

p a,N.S. 

atm 

p imp 

atm 

U a 

m/sec 

U D 

m/sec 

p a,D 

P b 

M = — 
D % 

10 

313.16 

1632 

1077 

8.415 

14.87 

12.60 

751 

1277 

1.700 

3-59 

20 

313.16 

2524 

1485 

12.96 

22.65 

21.25 

887 

1575 

1.777 

4.42 

29.6 

313.16 

2864 

1625 

14.72 

25.34 

24.89 

919 

1668 

1.8l4 

4.68 

40 

313.16 

2919 

1650 

14.96 

25.82 

25.26 

931 

1686 

1.811 

4.72 

50 

313.16 

2717 

1573 

13-68 

24.34 

22.31 

932 

i64o 

1.760 

4.58 


TABLE 3- 


PROPERTIES OF A CARBON MONOXIDE - AIR DETONATION WAVE, 


p =5 atm 
b 


Molar 
per cent 
fuel in 
mixture 

T b 

°K 

T 

a,D 

°K 

T 

a,N.S. 

°K 

P a,D 

atm 

P a,N.S. 

atm 

P. 

imp 

atm 

u 

a 

m/sec 

U D 

m/sec 

p a,D 

P b 

11 

.P 

10 

313.16 

1632 

1077 

42.07 

74.33 

63.00 

751 

1277 

1.699 

3.59 

20 

313.16 

2563 

1501 

65.89 

114.8 

108 .6 

889 

1586 

1.784 

4.45 

29.6 

313.16 

2982 

1673 

76.01 

131.3 

128.5 

936 

1698 

1.813 

4.76 

4o 

313.16 

3008 

1688 

75.82 

132.7 

127.1 

953 

1709 

1.794 

4.79 

50 

313.16 

2728 

1579 

68.30 

122.3 

110.9 

939 

1644 

1.751 

4.59 





TABLE 4 . PROPERTIES OF A CARBON MONOXIDE - AIR DETONATION WAVE, p = 25 atm 

b 


Molar 
per cent 
fuel in 
mixture 

T b 

°K 

T a,D 

°K 

T 

a,N.S. 

°K 

p a,D 

atm 

P a,N.S. 

atm 

P* 

imp 

atm 

U a 

m/sec 

U D 

m/sec 

p a,D 

P b 

*7* 

df 

11 

10 

313.16 

1632 

1077 

210.4 

371.7 

315.0 

751 

1277 

1.699 

3.59 

20 

313.16 

2583 

1510 

329.9 

578.6 

541.1 

897 

1592 

1.775 

4.47 

29.6 

313.16 

3090 

1717 

389.5 

677.3 

658.4 

955 

1724 

1.806 

4.83 

40 

313.16 

3066 

1714 

381.4 

675.9 

633.6 

971 

1725 

1.777 

4.83 

50 

313.16 

2732 

1581 

341.4 

612.5 

553.4 

941 

1645 

1.749 

4.60 


TABLE 


PROPERTIES OF A CARBON MONOXIDE - AIR DETONATION WAVE, 


p = 100 atm 
b 


Molar 
per cent 
fuel in 
mixture 

°K 

T a,D 

°K 

T 

a,N.S. 

°K 

P a,D 

atm 

P a,N.S. 

atm 

P. 

imp 

atm 

U 

a 

m/sec 

U D 

m/sec 

p a,D 

P b 

U 

\-T 

b 

10 

313.16 

1632 

1077 

841.4 

1487 

1260 

751 

1277 

1.699 

3.59 

20 

313.16 

2589 

1514 

1312 

2322 

2136 

905 

1595 

1.762 

4.48 

29.6 

313.16 

3173 

1750 

1583 

2773 

2666 

971 

1744 

1.797 

4.89 

1(0 

313.16 

3093 

1725 

1534 

2725 

2545 

976 

1732 

1.774 

4.85 

50 

313.16 

2733 

1582 

1365 

2451 

2209 

1 

942 

1646 

1.747 

4.6o 





TABLE b 

EQUILIBRIUM CONSTANTS 


T 

°K 

Atm. 

K. 

°2 

Atm. 

So 

Atm. 


-6 

-6 

-3 

1500 

4.843 x 10 

4.628 x 10 

3.O87 x 10 


-3 

, -4 

-2 

2000 

1.310 x 10 

7.37^ x 10 

1.920 x 10 


-2 

-2 

, , -2 

2500 

3.627 x 10 

1.564 x 10 

5.744 x 10 


, -1 

-l 

-1 

3000 

3.242 x 10 

1.205 x 10 

1.190 x 10 



„ -1 

-1 

3500 

1.517 

5.198 X 10 

1.997 x 10 




, -1 

4000 

4.759 

1.559 

2.941 x 10 




-1 

4500 

11.43 

3 .664 

3.957 x 10 




, n - 1 

5000 

22.77 

7.271 

4.998 x 10 


66,7^7 cal/mole CO 
58,532 cal/mole 0 2 
21,600 cal/mole NO 









TABLE 7 


T 


HEATS OF FORMATION 


A H„ 




T 

°K 

co 2 

CO 

°2 

0 

1 " 

n 2 

NO 

298.16 

-158.7954 

-44.6155 

0 

99.8153 

0 

36.6859 

300 

-157.7929 

-44.3197 

0.02266 

99 .2197 

0.02216 

36.4839 

400 

-117.1620 

-32.3613 

0.91113 

75 .0662 

0.89401 

28.2603 

500 

- 92.69516 

-25.1785 

1.46526 

60.5665 

1.42236 

23-3355 

600 

- 76.32368 

-20.3781 

1 . 854 it 8 

50.8971 

1.78348 

20.0650 

700 

- 64.58725 

-16.9379 

2.14861 

43.9881 

2.05086 

17.7426 

1000 

- 43.32795 

-10.6938 

2.73212 

31.5483 

2.58177 

13.6129 

1500 

- 26.61104 

- 5.75287 

3.25747 

21.8691 

3.08103 

10.4794 

2000 

- 18.16171 

- 3.23862 

3.56067 

17.0279 

3.37749 

8.95096 

2500 

- 13.05162 

- 1.71094 

3.76959 

14.1237 

3.57631 

8.05059 

3000 

- 9.62325 

- 0.683082 

3.93051 

12.1885 

3.7199^ 

7.45841 

3500 

- 7.16045 

- 0.056999 

4.06172 

10.8084 

3.82800 

7 . 0409 o 

4000 

- 5.30190 

0.615081 

4.17082 

9.77631 

3.91357 

6.73253 

4500 

- 3 .84864 

1.05183 

4 . 264 l 8 

8.97617 

3.98225 

6.49661 

5000 

- 2.68039 

1.40334 

4.34421 

8.33927 

4.03972 

6.31331 
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FIG. 2 PRESSURE RATIOS ACROSS CARBON MONOXIDE - AIR DETONATION WAVES 
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FIG. 3 GAS TEMPERATURES BEHIND CARBON MONOXIDE - AIR DETONATION WAVES 
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FIG. 7 HOLE FRACTIOUS OF SPECIES; INITIAL PRESSURE = 1 ATMOSPHERE 
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9 MOLE FRACTIONS OF SPECIES; INITIAL PRESSURE = 25 ATMOSPHERES 
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